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Abstract 

Nonlinear random vibration under excitations of both Gaussian and Poisson white noises 
is considered. The model is based on stochastic differential equations, and the corresponding 
stochastic integrals are defined in such a way that the energy conservation law is satisfied. It is 
shown that Stratonovich integral and Di Paola-Falsone integral should be used for excitations of 
Gaussian and Poisson white noises, respectively, in order for the model to satisfy the underlining 
physical laws (e.g., energy conservation). Numerical examples are presented to illustrate the 
theoretical results. 

Keywords: Random vibration, nonlinear systems, Poisson noise, Gaussian noise, Stochastic 
differential equations, stochastic integrals. 

1 Introduction 

Differential equations are extensively used in modeling dynamical systems in science and engi- 
neering. When dynamical systems are under random influences, stochastic differential equations 
(SDEs) may be more appropriate for modeling. The solutions of SDEs are interpreted in terms of 
stochastic integrals [91 [6] . 

Dynamical systems subject to Gaussian white noise are often modeled by SDEs with Brownian 
motion, and the solutions are in terms of the Ito integral |1CH El [6]. Although the Ito integral is 
self-consistent mathematically, it is not the only type of stochastic integrals that can be constructed 
to interpret an SDE. Other stochastic integrals, such as the Stratonovich integral [10\ [9| [6], have 
also been used to interpret an SDE as a stochastic integral equation. There is no right or wrong 
choice when choosing either Ito or Stratonovich integrals in interpreting SDEs mathematically, 
since the two integrals are equivalent and can be converted into each other, provided that the 
integrand satisfies certain smoothness conditions |1CH El [6]. However, these stochastic integrals 
have different definitions, and one may be more directly related to a practical situation than the 
other. While Ito integral is a reasonable choice in many applications including finance and biology 
[9], Stratonovich integral is believed to be more appropriate in physical and engineering applications 
[lOj . Stratonovich integral has an extra term comparing with the corresponding Ito integral: the so- 
called correction term \10\ [T2] . Some authors [13] attribute this correction term to the conversion 
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from physical white noise to ideal white noise. This explanation is not necessarily convincing 

nam. 

Dynamical systems driven by non-Gaussian white noise, especially Poisson white noise, have 
attracted a lot of attention recently. Correction terms for converting Ito SDEs to Stratonovich 
SDEs with Poisson white noise are presented in [IJ [2] . Although these correction terms have been 
accepted widely, there are some confusions [3l 0] . 

In this paper, we consider nonlinear random vibration under excitations of either Gaussian or 
Poisson white noises, modeled by appropriate stochastic differential equations. The main objective 
of this paper is to explain the correction terms in both Gaussian and Poisson white noise cases, 
from a physical perspective. We will show that the correction terms are natural consequences 
of fundamental physical laws satisfied by the vibration system. Note that conventional spectral 
analysis methods [8], which have found extensive applications in random vibration analysis, are 
not applicable in this case due to the nonlinearity of the system. 

To this end, we consider a vibration system as a mass-spring-damping oscillator with random 
excitation 

mx(t) + kx(t) = g(x(t),x(t)) + f(x(t),x(t))L(t), (1) 

where m represents the mass, k is the stiffness coefficient of the spring, x(t) is the displacement 
depending on time t, and x(t) is the velocity. g(x(t),x(t)) and f(x(t),x(t))L(t) represent the 
generalized force terms, which may originate from external or parametric excitations. L(t) is a 
noise term defined as the formal derivative of some stochastic process 

L(t) = bB(t) + cC(t), (2) 

where b and c are constants, Bit) is a Gaussian process, and C(t) is some compound Poisson 
process, which is expressed as 

N(t) 

C(t) = ]T RiU(t - U). (3) 
i=i 

In Eq. ([3]), N(t) is a Poisson process with intensity parameter A, U(t — ti) is a unit step function 
(a Heaviside function) at ti, Ri is a random variable representing the z-th impulse. It follows from 
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© that 

L(t) = bB(t) + cC(t), (4) 
where B(t) is the Gaussian white noise, and C(t) is the Poisson white noise expressed as 

N(t) 

C{t) = Y J ^-U). (5) 

i=i 

Note that ([!]) expresses a general noise model including the Gaussian white noise (6 ^ 0,c = 0), 
the Poisson white noise (6 = 0,c ^ 0), and the combined Gaussian and Poisson white noise (b ^ 
and c ^ 0). 

The second-order equation ([1]) can be rewritten as a system of SDEs 

M)\Jo i)(*m) it+ L( » ) it+ L( ° | dL(t , (6) 

Since L(i) is non-differentiable almost everywhere, (jSj) cannot be interpreted in the framework of 
classical calculus. Thus the solution of © is interpreted with a stochastic integral, 

x{>) ) = h) + r ( ° k 'V-wU+irf i ds 

£(t)/ \±(0)/ -/o Oj y±( s )J m h \g(x(s),x(s)) 



1 ft / 



■/?? 



o \/(x(s),x(s)) 



dL( S ). (7) 



Defining y(t) = I I , A = l I and using the variation of parameters formula, the 
solution to Eq. ([!]) can also be rewritten as [9], 



^ ds+- f e A{t ~ s) I I dL(s), (8) 



g(x(s),x(s))J m Jo \f(x(s),±(s)) 

where yo = = ^ s the initial condition. It can be shown that ([7]) and flS]) are 



±(0) I Vxq 



equivalent [9]. 
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It is straightforward to verify that 



, At AH 2 

I+ T! + ^r + 



-wsin(wi) cos(wi) 



(9) 



where 



— . Substituting ([9]) into flSj), we get 



. . sin(wt) 

= cos(wi)2;o H xq + 



+ 



Ul 

sina;(t — s) 



sin u>(t 



moo 



-g(x(s),x(s))ds 



x{t) 



+ 



raw 

uj sin(wt)xo + cos(wt)i;o + 
* coso;(t — s" 1 



f(x(s),x(s))dL(s), 

X* cosu;(t — s] 
m 



(10) 



g(x(s),x(s))ds 



m 



f(x(s),x(s))dL(s). 



Note that the stochastic integrals in Eqs. (|T|) and (]10p are yet to be defined. As stated earlier, 
the stochastic integrals which can be used to interpreted SDEs may not be unique. The question 
is that which stochastic integral will lead to the solution that is consistent with the physics of 
the system. One possible answer is to compare solutions to the SDEs with the corresponding 
experimental results. However, this method may be impractical in many cases due to the high cost 
of performing the experiments, as one needs highly accurate data from sufficiently large number of 
samples in order to resolve the subtle difference in the theory. In this paper, to construct a SDE 
model that is physically relevant to the real system, we propose to apply a stochastic integral such 
that the fundamental physical law (e.g., energy conservation) is satisfied. 

This paper is organized as follows. In Sec. [21 starting from the energy conservation law, we define 
the stochastic integral that is suitable for the SDE model of the nonlinear random oscillators. The 
relationship between the proposed models and the existing models is discussed in Sec. El Numerical 
methods with an illustrative example are presented in Sec. [H 



2 Stochastic integrals for nonlinear oscillators under noise excita- 
tion 

Since there are multiple forms of the stochastic integrals that can be constructed from the SDE, 
we define the stochastic integral such that the fundamental physical laws are satisfied. As for the 
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nonlinear oscillators described by the SDE ([T]) , we expect the energy- work conservation be satisfied 

"1 .o,., 1 



2 mx 2 (t) + -kx 2 {t) 



^mx 2 (0) + ^kx 2 (0) 



f 

Jo 



g(x(s),x(s)) + f(x(s),x(s))L(s) dx(s), 



where ^mx 2 (t) + ^kx 2 (t) represents the total mechanical energy of the system at time t, and the 
integrand in the right hand side is the forcing term of (pQ) . Equation (jlip expresses that the change 
in the total mechanical energy is equal to the work done by the external forces. 
Writing (llip in the form of stochastic integral, we have 



Knx 2 {t) + hx 2 {t) 



^mx 2 (0) + ^kx 2 (0) 



t rt 

g(x(s),x(s))x(s) ds + / f(x(s), x(s))L(s)x(s) ds 
Jo 

t ,t 

g(x(s), x(s))x(s) ds + / f(x(s), x(s))x(s) dL(s). 
o Jo 



(12) 



As stated before, the stochastic integral with respect to L(t) should be defined such that the 
solution of ([T]) satisfies the energy conservation law (|12p . It follows from ^ that a stochastic 
integral with respect to L(t) can be decomposed into two terms: stochastic integral with respect 
to B(t) and stochastic integral with respect to C(t). We define the two terms in the next two 
subsections. 



2.1 For Gaussian white noises 

Assume 6=1 and c = 0, then it follows from ([2]) that the stochastic process Lit) reduces to a 
Brownian motion, and (HOD and (1121) become 



x(t) 



x(t) 



sm(ujt) 
cos{ujt)XQ H xq + 



I' 



sin u>(t 



moo 



sinwft — s) „ 

i if (x s), x s dB(s), 

moo 



-g(x(s),x(s))ds 



+ 



Jo S ^ 

■ r \ i n • f ' cosu;(t — s 

-u) smlutjxo + cos(u)t)XQ + / 

X' 



(13) 



g(x(s),x(s))ds 



COS L0(t 



-f(x(s),x(s))dB(s). 
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and 



\mx 2 {t) + \kx 2 {t) 



\mx 2 (<d) + \kx 2 {0) 



t f t 
g(x(s), x(s))x(s) ds + / f(x(s), x(s))x(s) dB(s), 
Jo 



(14) 



respectively. 

There are two types of stochastic integral extensively used for SDEs driven by Brownian motions: 
Ito integral and Stratonovich integral. Throughout this paper, we use to denote Ito calculus, 
and V for Stratonovich calculus. In the sense of Ito, (fl~3l) and (fTH can be written as 



and 



. sin(tjt) . r t smoj(t — s) 

x(t) = cos{ujt)XQ H xo + / g(x{s), x{s)j ds 

uj Jo 

^ sina;(t — s) 



+ 



moo 



-f(x(s),x(s))*dB(s), 



x{t) 



r* coswft — s 

—uj s\n(ijjt)xQ + cosiijjtjXQ + I 

Jo m 

+ r COSU{t ~ S) f(x(s),x(s))*dB(s). 
Jo m 



g(x(s),x(s))ds 



(15) 



^mx 2 (t) + ^A;x 2 (t) 



^mx 2 (0) + ^kx 2 (0) 



t rt 

g(x(s), x(s))x(s) ds + / f(x(s), x(s))x(s) * d.B(s), 
o j o 



(16) 



respectively. In the sense of Stratonovich, (|13p and (|14p can be written as 



and 



. . sin wt . f* smw t-s 

= cos{ujt)XQ H xo + / g(x{s), x{s)) ds 

UJ Jo 
J^ 1 sina;(t — s 



x(t) 



+ 



muj 



muj 

'-f(x(s),x(s))odB(s), 



... P cosu(t — s) 

-uj smlujtjxo + cos(ujt)XQ + / 

■a) m 



g(x(s),x(s))ds 



cos uj(t — s) 



m 



f(x(s),x(s))odB(s). 



(17) 



^mx 2 (t) + ^kx 2 {t) 



imi 2 (0) + ^fcc 2 (0) 



i rt 

g(x(s), x(s))x(s) ds + / f(x(s), x(s))x(s) o dB(s), 
Jo 



(18) 
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respectively. Provided that the function / is sufficient smooth, the solutions in Stratonovich inte- 
grals, (|17p and (|18p . can be converted into the following forms with Ito integrals [9] 

X t sinw(i — s) 



xyt) = cosyuitjxQ H xq 



g(x{s),x{s)) + —f(x(s),x(s))U[x(s),x(s)) 



els 



sinw(< — s) 



f(x{s),x(s))*dB(a), 



r 

x(t) = — ui sm(uit)xo + cos(uit)io + I 

X 1 cos ui(t — s) 
m 



cos uj(t — s) 



m 



g(x(s),x(s)) + —f(x(s),x(s))f ± (x(s),x(s)) 



ds 



f(x(s),x(s))*dB(s), 



(19) 



and 



^m± 2 (t) + ^kx 2 (t) 



g(x(s),x(s)) + — f(x(s),x(s))fa;(x(s),x(s)) ) x(s) + ^—f 2 {x(s),x(s)) 
2m J 2m 



+ / f(x(s),x(s))x(s)*dB(s). 

'0 



ds 



(20) 



As shown in the Appendix, the solution (|19p satisfies the energy- work relation (|20p . suggesting 
that when the randomness is modeled in sense of Stratonovich, the energy-work conservation law 
is satisfied. On the other hand, in a similar procedure as in the Appendix, it can be shown that the 
energy- work law (|16p contradicts with the solution (|15p . Therefore, Stratonovich integral instead 
of Ito integral should be used so that this nonlinear random oscillator model satisfies the energy 
conservation law. 

This implies that when Gaussian noise is present in this nonlinear vibration system, the SDE 
model should be interpreted in the sense of Stratonovich stochastic integral, but not in the sense 
of Ito stochastic integral. 



2.2 For Poisson white noises 

When 6 = and c = 1, the stochastic process as expressed in ([2]) reduces to a compound Poisson 
process. Note that the jump size of C(s) at time s can be expressed as AC(s) = C(s) — C(s-), 
where C(s-) is the left limit of C(s) at s. Suppose C(s) jumps at times U (i = 1, 2, • • • ), then the 
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solution ([7]) can be written as 
x(t) = x(0) + x(s)ds, 

k p* i p* i N ® c u 

(21) 

x(t) = x(0) I x(s)ds + — I g(x(s),x(s))ds + — V / /(x(s),x(a)) dC(s), 

where N(t), as shown in ([3]), represents the number of jumps upto time t. 

In the following, we shall derive the stochastic integral with respect to jumps such that the 
energy conservation law is satisfied. First, let's examine the changes in the system at i-th jump 
occured at time ti (1 < % < N(t)). From (|2ip . the displacement x is continuous while the velocity 
x undergoes an jump given by 

x(ti) = x(ti-), 

! PU ( 22 ) 

xiti) = x{ti-) + — I f(x(s),x(s))dC(s). 
< m Jti- 

The change in the total energy (|12|) due to the i-th jump is that in the kinetic energy given by 

lrnx 2 (ti) - imi 2 (4 r ) = f ' f(x(s), x(s))x(s) dC(s), (23) 

J ti — 

due to the continuity of the displacement x across an jump. 

If the integrals with respect to jumps are defined in sense of Ito, then (|22p and (|23p becomes 



x(U) = x(U~) + -f(x(ti-),x(ti-))AC(ti)- (24) 
m 

and 

^m± 2 (ti) - ^m± 2 {ti-) = f(x(U-),±(U-))x{U-)AC(U), (25) 

respectively. Since AC(ij) ^ 0, it is obvious that (|24p contradict with (|25p . which indicates that the 
energy conservation law cannot be satisfied when the integrals with respect to jumps are interpreted 
in sense of Ito. In the following, we shall show that the integrals should be interpreted as some 
kind of Riemann integral on the imaginary path along the jump to satisfy the energy conservation 
law. 

Let x(ti, r) be the value of x(s) at time ti if C(s) jumped from C(tj— ) to r. Then x(ti, C(ti— )) = 
x(ti— ) and x(ti,C(ti)) = x{tj). With the integrals being interpreted as the Riemann integral on 
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the imaginary path along the jump, the energy-work law (|23p can be written as 

-mx 2 (t i ,C(ti))--mx 2 {t i ,C(t i -))= / f(x(U), r))x(U,r) dr. (26) 
2 2 Jc{u-) 

and the solution (f22|) becomes 

1 fC(ti) 

x(t i ,C(t i ))-x(U,C(U-)) = - f(x(U),i(U,r))dr. (27) 

m JC(ti-) 

Since the jump size can be any value, it follows from (j26|) and (f27|) that for any A 6 R, it is true 
that 

\mx?{t i ,\)-\m£{t il C{t l -))= [ f(x(t i ),x(t i ,r))x(t i ,r)dr, (28) 

and 

x(ti,\)-x(ti,C(ti-)) = — [ f(x(ti),x(ti,r))dr. (29) 

m yc(ti-) 

Taking derivatives of both sides of (I28p and (I29D with respect to A, respectively, we get the identical 
ordinary differential equation(ODE) 

^rx(ti, A) = -f(x(U), x(U, A)). (30) 
dA m 

Therefore, the energy conservation law is satisfied. 

Using the fact that x(ti,C(ti—)) = x{ti~) and x(ti,C(ti)) = x(ti) , it follows from (f30j) that 

x(ti) = x{U-) + -Yi(AC(ti)), (31) 
m 

where li(AC(tj)) is determined by the initial or terminal value problem of the ODE 

' ^-Y l (\) = f(x(U),Y i (\) + ±(t i -)), 
dA (32) 

Yi(0) = 0. 

Note that in fl32]), < A < AC(U) for AC(U) > or AC(U) < A < for AC(tj) < 0. Comparing 
the original solution expression (|22p with the new formula (|3ip . it can be seen that the last term 
in (f22l) should be defined as Yj(AC(tj)), and hence (f2"T|) should be interpreted as 

x(t) = x(0) + J * x(s) ds, 

i ft i iV(t) (33) 

x{t)=x(0) + - I g(x(s),x(s))ds + - Yi(AC(ti)), 
m Jo m i = \ 
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where Y;(AC(^)) is the solution to the ODE (|32j). 

As it will become clear in Sec.O this implies that when pure jump noise, such as Poisson white 
noise, is present in this nonlinear vibration system, the SDE model should be interpreted in the 
sense of Di Paola-Falsone stochastic integral [H [2] . 

2.3 For combined Gaussian and Poisson white noises 

When both 6^0 and c / 0, the excitation is a combined Gaussian and Poisson white noise. 
Combining the results in the subsections 12.11 and 12.21 we find that, in order to satisfy the energy- 
work conservtion law, one has to interpret the stochastic integrals with respect to Brownian motions 
as Stratonovich integrals, and the integrals with respect to jumps as DiPaola-Falsone integrals. 
Therefore, the solution to (pQ) is given by the expression Q, where the stochastic integral is defined 
as 



where AL(ij) = c[C(ti) — C(t{—)] = cAC(t{). Recall that, in (|34l) . V denote integrals in the 
Stratonovich sense, N(t) represents the number of jumps up to time t, and Yj(AL(ij)) is the 
solution to the ODE (|3"5j) . where A takes value of < A < AL(U) for AL(U) > or AL(U) < A < 
for AL(ti) < 0. 

3 Relationship with the existing models 

In this section, we shall show that the correction term Yj(AL(ij)), as given by the solution to the 
ODE (|32p . is consistent with the one proposed in the work by Di Paola and Falsone [H[2]. 

When / is Lipschitz continuous, it is easy to check that the solution of (I32p exists and is unique. 
If / is assumed to be smooth, then Yj(AL(tj)) is analytic with respect to AL(ij). 

Using Taylor expansion, 




t ft Jv w 

f(x(s),±(s)) dL(s) = b / f(x(s),±(s)) odB s + ^ Fj(AL(ij)), 



(34) 



Yi(AL(ti)) = Yi(0) + — y,(A)| A=0 (AL(^)) + 



l__d^ 
2!dA 2 



y J (A)| A=0 (AL(t i )) 2 + ... . 



(35) 



It follows from (|32p that for any n > 1 




(36) 
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where 



Substituting ([55]) into (f3"5j) . and using the fact that 1^(0) = 0, we get 

Y t (AL(U)) = f; / 0) (*(*t),*(*t-)) {Am) J (37) 

V^C*),*^-)) = /(sfej.ife-)) 
/<%(*),*&-)) = g/ M) (g(*0,^-) + A) fori > 2 . 

OA A=0 

Thus the correction term given by (|37p is exactly the same as the one proposed in [TJ [2] . 

We have shown that the correction term li(AL(tj)) can be obtained in two ways: solving the 
initial value problem to the ODE (|32p or computing the expansion (|37|) . Note that the former 
approach of solving ()32|) is applicable under much more general condition than the latter one of 
evaluating the infinite series (|37p . because the existence of the solution to the ODE only requires 
f(x,y) to be integrable but (f3T|) demands f(x,y) to be infinitely differentiable. 



4 Simulation examples 

Solutions of the SDE for a nonlinear oscillator ([1]) , defined by (J7J) and ()34|) , can hardly be obtained 
with analytical methods. In this section, the SDE ([1]) is numerically solved to verify the conclusion 
obtained in section 2. Consider the case with both Gaussian and Poisson white noises, dL = 
bdB + cdC, with the compound Poisson process given by ([5]). 

The numerical procedure of the SDE for a nonlinear oscillator (JTJ) defined by ([7]) and (|34p is as 
follows. On each time subinterval U-± < t <U, i.e. when no jumps occur, ([T]) becomes 

dx(t) = x(t), dt 

k 1 b (39) 

dx(t) = x(t) dt H g{x(t),x(t)) dt H f(x{t),x(t)) o dB. 

mm m 

The above equation can be converted into Ito SDE and then computed by conventional algorithms 
for Ito SDEs, such as Euler method, Milstein method, or other algorithm of high-order accuracy 
based on stochastic Taylor expansion [?]• At the time ti when an jump occurs, it follows from ([7]) 
and (El) that 

x(ti) = x(U-), 

(40) 

x(U) = x(U-) + ±Yi(AL(U)), 
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where K;(AL(ij)) is obtained by solving the deterministic ODE (|32|) using Runge-Kutta or multistep 
methods. 

Consider the following stochastic Duffing-van der Pol equation 

x + (1 + x 2 )x + x + x 3 = ±L(t), (41) 

with the initial condition x(0) = 2 and x(0) = 0. The SDE (|4ip can be written in form of (pQ) with 
m = 1, k = 1, g{x,x) = —(1 + x 2 )x — x 3 and f(x,x) = x. In the simulation, we take L{t) as in 
([2]) with b = c = 1, i.e. L(t) = B(t) + C(t), where C(t) is a pure jump process given by ([3]) with 
N(t) being a Poisson process with intensity parameter as A = 3.4 and Ri {i = 1, ■ ■ ■ ,N(t)) being 
random numbers of the standard normal distribution. 

Case 1 

In this case, (f4T|) or ([1]) is interpreted by ([7]) and (i34"|) . In the simulation, to integrate ([7]) and 
(I34p numerically, we use Euler's method to advance fj39|) when no jumps occur, while evaluate (|40j) 
by solving (|32p with Euler's method when jumps arrive. Note that to apply Euler's methods, the 
stochastic integral in (|39p need to convert into Ito integral. The step size of Euler's method for 
solving both (|39h and (|32|) is At = 0.0001. Figure Q] shows a sample path of the driven process 
L(t), and Figs. [2] and [3] show the numerical solution of the displacement x(t) and the velocity x(t) 
respectively, corresponding to the path shown in Fig. [TJ Based on the numerical solutions shown 
in Figs. [2] and [3l we compare in Fig. [4] the energy increment AE(i), which is defined by 

AE = ^ [x 2 (t) + x 2 (t) - x 2 (0) - ± 2 (0)] , (42) 

with the work done WK(i) defined by / g(x, x)x(s) ds + f(x, x)x(s) dL(s) where the stochastic 

Jo " Jo 



integral is taken in the sense similar to (|34|) . i.e. 



WK(t)= / g(x,x)x(s)ds+ / f(x,x)x(s)dL(s) 
Jo Jo 



/ g(x,x)x(s)ds+ / f(x,x)x(s)o dB s + Vl^AL^)), (43) 
Jo Jo ~i 



1=1 



where AL(tj) = c[C(ti) — C(tj— )] = cAC(i,), ! o' denote integrals in the Stratonovich sense, N(t) 
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is the number of jumps upto time t, and YJ(AL(ij)) is the solution to the following ODE, 

' ^YiW = f(x(ti)M^+i(U-))(YiW+i(t i -)), 
Yi(0) = 0. 

with A taking value of < A < AL(^) for AL(t;) > or AL(ij) < A < for AL(U) < 0. It can be 
seen clearly from Fig. [J] that the energy increment of the system agrees with the work very well, 
indicating that the energy conservation law is satisfied. 

Case 2 

In this case, (|41[) (or ([1]) is interpreted by using Ito stochastic integrals. Now (j39|) and (|40p 
become 

x(t), dt 

k 1 b (45) 

x(t) dt H g(x(t),x(t)) dt H f(x{t),x(t)) ★ dB, 

mm m 

and 

x(ti) = x(ti-), 

(46) 

x(ti) = x{ti-) + ±f(x(s-),x(s-))AL(U), 

respectively. In the simulation, the driving process L{t) and all the simulation parameters are 
taken the same as in case 1. Figures Oand O present the corresponding numerical solutions of the 
displacement and velocity, respectively. Comparison of the energy increment, defined by (|42p . and 
the work done, now defined by 

,* ,t m 

WK(t)= g(x,x)x(s)ds + b f(x, x)x(s) * dB s + V f(x(s-), i(s-))i(s-)AI(t,), (47) 
Jo Jo i=1 

is presented in Fig. We can see clearly from Fig. [7] that there is significant difference between the 
energy increment and the work done by the force. Note that all the curves in Fig. [7] tend to have 
a very small variance in the time span 0.84 < t < 1. This is the consequence of the fact that the 
velocity is very small for 0.84 < t < 1, as shown in Fig. [6l Since f(x(t),x(t)) = g(x(t),x(t)) = x(t), 
it follows from (|42p . (|45p and (|47p that both the energy increment and the work done change slowly 
for very small velocity x(t). 
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By comparing Fig. [7] with Fig. we can see that Stratonovich integral and Di Paola-Falsone 
integral should be used for excitations of Gaussian and Poisson white noises, respectively, in order 
for the model to satisfy the underlining physical laws. 



A Appendix: Proof of the energy-work law ( 1201) from the solution 

<m 



For cosmetic purpose, we introduce the following simplified notations: g(s) = — g(x(s), x(s)), 

m 

f(s) = — f(x(s),x(s)), fx(s) = — fx(s)( x ( s ), x( s ))- Moreover, in this Appendix, all the stochastic 
m m w 

integrals with respect to Brownian motions are in sense of Ito (we have dropped * notation). Then 
the energy-work law (f20|) is equivalent to 



i[x 2 (t) +W V(t)]-I[i 2 (0)+a;V(0)] 



o 



1 



9(s) + -f(s)f x (s)jx(s) + -f(s) 



1 



ds+ [ f(s)x(s)dB(s). 
Jo 



(48) 



Next, we show the solution given in (|19p satisfies the energy- work law (|48p . 

Denote the right and left hand sides of (j48l) as RHS and LHS, respectively. Substitute ([19 
into the left hand side of (08]), we get 

2 



g(s) + lf(s)h( s ) 



cos u>(t — s) 



t \ 2 

sin(w (t - s))/(s) dB(s) 



ds) + 



1 

2 vo 

t x 2 

cos(w(t- s))f(s)dB(s 



+ (cj cos(uit)xo + sin(wi)±o) / sin uj(t — s) 



ds 



+ (lu cos(wi)xo + sin(wt)io) / sm(w(t — s))f(s) dB(s) 

Jo 



+ ( J sinuj(t — s) 



sm(u(t - s))/(s) dB(s) 



+ (— ui sm(uit)xo + cos(u>t)io) / cosw(t — s) 



g(s) + 5 /(*)/x(«) 



ds 



+ (— oj sm(wt)x + cos(cot)x ) / cos(uj(t — s))f(s) dB(s) 

Jo 



+ ( J cos(u(t — s) 



l(s) + -f(s)U(s) 



(Is 



cos(w(i - s))/(s) dB(s) 



(49) 
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Substituting (fTUj) into the right-hand side of (|48p. we get 



t 

RHS = I (g(s) + f(s)fx(s)) (—uj sin(u!s)xQ + cos(cjs)xo) ds 
'o 



t rs 



+ 



+ 




'0 JO 

r-t r s 



cos(w(s -p)) g(s) + 



'o J o 

ft 



gip) + 



f(p)f±(p) 



dpds 



d£(p)d s + iy* f(s)ds 



+ / /(s) (— w sin(a;s)a;o + cos(o;s)i;o) dB(s) 
J o 

+ jf jT cos(a;( S - p))/( S ) ^(p) + IM£M \ dpdB{s) 

+ I f co S (co(s-p))f(s)f(p)dB(p)dB(s). 
Jo Jo 

To prove LHS in (|49p is equal to RHS in (|50p . we claim the following facts 



2 



1 
2 

+ 



^ sin(u,(t - s)) (g(s) + MM^j ds 
\os(u(t-s))(g(s) + f{s)f 2 b{s) ) ds 



/ [ S cos(u(s-p))f(s)f(p)dB(p)dB(s) + \ [ f 2 (s)ds 
Jo Jo z Jo 

= X - Qf * sm(u(t - s)f(s) dB(s)j + i Qf cos( W (t - d£(s)) 



cos{uj{s-p) g{s) + 




o Jo 



+ 



f(p)dB(p)ds 



t rs 



JO 



cos(u}(s-p) gip) + 



f(p)fx(p) 



f(s)dpdB(s) 



^ cos(u)(t - s)) yg(s) + 
t 



f(s)h(s) 



ds I cos(w(f - s))f(s) dB(s) 



t 



- f sm(u{t - s)) (g(s) + MM*T\ ds J sin(o;(< - s))f(s) dB(s) 
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g{s) + ^ - ) (— oj s'm(ujs)xQ + cos(a;s)xo) ds 







(— uj s\n{ut)xQ + cos(wt)i;o) J cos(uj(t — s)) (g(s) + ^ s ^ x ^ ) 

l*~b / £ t \ £ ( \ 

+ (a; cos{ojt)xQ + sin(a;t)xo) / sin(a;(t — s)) I g(s) -\ — ) ds, (54) 



and 



o 



o 



t 

f(s) (— uj sin(ws)xo + cos(a;s)i;o) dB(s) 

ct 



2 



= (— uj sin(cjt)xo + cos{ujt)xo) / cos(w(i — s))f(s) dB(s) 

Jo 

+ (uj cos(ujt)xQ + sin(u;t)xo) / sin(a;(i — s))/(s) cLB(s). (55) 

J o 

One can easily see that ([54"|) and ([55]) are true by using the trignometric identities 
cos(ws) = cos(wi) cos(w(t — s)) + sin(cjt) sin(u;(i — s)), 

and 

sin(ws) = sin(o;i) cos(a;(i — s)) — cos(wt) sin(w(t — s)), 



In the following, we give the proofs for (I5ip and (I52p . The proof of (I53D is similar to those for ([51 
and (152(1 and is not given here. 

To prove (j5ip is true, we rewrite the right-hand side of (|5ip as double integrals 



* - -)) f + d*V + f /' c OS(W (t - .)) f ff ( fl ) + zifiMf) ) d ; ' 



o \ * / / vo 

t /.i 



sin( W (t - p)) sin( W (t - «)) ( ff (p) + f g(s) + ) dpds 



JO 

* ^COBWt-p)) C08(a,(t-*)) Up)+ fiP)f 2 i(P) ^ (g( S )+ fiS) -!; {Sy ) dpds 

^ r \osHs- P) ) ( g(p)+ mm} ( 9{s)+ mm^ dPds 

2 f 008(0,(3 -p)) + (,(.) + dpd, 



JO 



2 / V 2 



Similarly, to prove (I52p . we rewrite the right-hand side (|52[) as 

t \ 2 / f* x 2 

sin(w(t-s))/(s)cLB(s)) + / cos(w(t - s))/(s) dS(s) 



/ \J0 



/ / cos(o;(*-p))/(*)/(p)dB(p)dB(*). (56) 
Jo Jo 
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The integral domain for the right-hand side of (|56|) is a square given by A = {(s,p) Is 6 [0,t],p e 
[0, i]}. Decompose the square into three parts: = {(s,p) < s < p < t}, A2 = {(s,p)\0 < p < 
s < t}, and ^3 = {(s, s) |0 < s < t}, then the right-hand side of (|56p becomes 

^os(w( S -p))/( S )/(p)d£(p)d£( S ) = // cos(w(a-p))/(a)/(p)dB(p)dB(a). (57) 



jo J J 

A1+A2+A3 

Note that 

cos(w( S - p))f(s)f(p) dB(p) dB(s) =H cos(u;( S - p))f(s)f(p) dB(p) dB(s) 
Ai JJa 2 

1 [ S cos(a;( S - p))f(s)f(p) <LB(p) <LB(s), (58) 
jo 

and 

ft 



^jf cos(u(s-p))f(s)f(p)dB(p)dB( 



(s)= / /^( 5 )d S . (59) 







It follows from (|5Sj) and K5SJ that is true. 

Add Eqs. ([^T]l to together, we get LiJ5 = RHS, and hence (J2QJ) is true. 
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Figure 2: The evolution of the displacement x(t) as the solution of Duffing-van der Pol equation 
(|41|) defined by Stratonovich integral and Di Paola-Falsone integral, where the driving process L(t) 
is given in Fig. [TJ 
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Figure 3: The evolution of the velocity x(t) as the solution of Dufnng-van der Pol equation (|4ip 
defined by Stratonovich integral and Di Paola-Falsone integral, where the driving process L(t) is 
given in Fig. [TJ 
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Figure 5: The evolution of the displacement x(t) as the solution of Dufhng-van der Pol equation 
(|4ip defined by Ito integral, where the driving process L(t) is given in Fig. [TJ 
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Figure 7: Comparison of the change in total energy, defined by (|42p . and the work done by the 
force, defined by ([¥?]) ■ 
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